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Abstract 

Neutron star radii are primarily determined by the pressure of isospin asymmetric 
matter which is proportional to the slope of the nuclear symmetry energy. Available 
terrestrial laboratory data on the isospin diffusion in heavy-ion reactions at inter- 
mediate energies constrain the slope of the symmetry energy. Using this constraint, 
we show that the radius (radiation radius) of a 1.4 solar mass (Mq) neutron star is 
between 11.5 (14.4) and 13.6 (16.3) km. 
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With central pressures of 10^^ dynes/cm^ and gravitational binding energies 
of 10^^ ergs, neutron stars are among the most exotic objects in the universe. 
Impressive progress has been made on the observable properties of neutron 
stars, such as masses, radii, spectra, and rotational properties [1,2,3,4,5,6,7,8]. 
However, a precise neutron star radius measurement still eludes us. 

The theoretical understanding of these observable properties demands an un- 
derstanding of the relevant nuclear physics. For recent reviews, see Refs. [9,10,11,12,13,14,15,16,17]. 
The global properties of neutron stars: masses, radii, and composition, are de- 
termined by the Equation of State (EOS) of neutron-rich nucleonic matter, 
thus neutron stars are ideal astrophysical laboratories for investigating the 
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EOS. The EOS can be separated into two contributions, the isospin symmetric 
part (the EOS of nuclear matter, -Enuc) and the isospin asymmetric part (the 
nuclear symmetry energy, -Esym)- This separation is manifest in the relation 
E{p, 6) = -Enuc(p) + S'^Esymip), where p is the baryon density, 6 = {pn - Pp)/p 
is the isospin asymmetry, and p„ and pp are the neutron and proton densities. 
While many neutron star properties depend on both parts of the equation of 
state, the radius is primarily determined by the slope of the symmetry energy, 
-^sym(p) [9,10,11,12]. Unfortunately, our knowledge about the density depen- 
dence of the nuclear symmetry energy has been rather poor. Predictions of 
the symmetry energy by nuclear many-body theories vary significantly [18]. 
Because of its importance for neutron star structure, determining the den- 
sity dependence of the symmetry energy has been a major goal of the in- 
termediate energy heavy-ion community. Although extracting the symmetry 
energy is difficult because of the complicated role of isospin in the reaction 
dynamics, several observable probes of the symmetry energy have been sug- 
gested [19,20,21,22] (see also Refs. [23,24,25] for reviews). 



Some significant progress has been made recently in determining the density 
dependence of -Esym(p) using: (i) isospin diffusion in heavy-ion reactions at 
intermediate energies as a probe of the Esym{p) around the saturation den- 
sity [26,27,28,29,30,31], (ii) flow in heavy-ion collisions at higher energies to 
constrain the equation of state of nuclear matter [24], and (iii) the sizes of 
neutron skins in heavy nuclei to constrain £'sym(p) at sub- saturation densi- 
ties [32,33,17,34]. 



The observational determination of a neutron star radius from the measured 
spectral fluxes relies on a numerical model of the neutron star atmosphere 
and uses the composition of the atmosphere, a measurement of the distance, 
the column density of x-ray absorbing material, and the surface gravitational 
redshift as inputs. Many of these quantities are difficult to measure, thus the 
paucity of radius measurements. Current estimates obtained from recent x-ray 
observations have given a wide range of results. 



In this Letter, we combine recently obtained isospin diffusion data, informa- 
tion from flow observables, studies on the neutron skin of ^°^Pb, and other 
information to constrain the radius of 1.4 Mq neutron stars. 



We use the EOS corresponding to the potential [35,36,37] 
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Here r = 1/2 ( — 1/2) for neutrons (protons) with r ^ t', a = 4/3, and f'Ti^.f)) 
is the phase space distribution function. The incompressibihty Kq of sym- 
metric nuclear matter at po is set to be 211 MeV. Eq. 1 is an extension of 
the potential from Weike et al. [38,39] from symmetric to asymmetric mat- 
ter. The parameter x was introduced to mimic various predictions on Egymip) 
by using different many-body theories and effective interactions. It is a con- 
venient way to parameterize the uncertainty in the magnitude and density 
dependence of the symmetry energy while keeping the isospin-symmetric part 
of the EOS unchanged. The functions Au{x) and Ai[x) depend on x according 
to Au{x) = -95.98 - x^,and Ai{x) = -120.57 -|- x-^ such that the same 
saturation properties of symmetric matter and a value of £^sym(po) = 32 MeV 
are obtained. 

The isoscalar potential estimated from {Uneutron + Uproton)/'^ agrees very well 
with predictions from variational many-body theory [40]. In addition to be- 
ing useful at the lower energies discussed here, the underlying EOS has been 
tested successfully against nuclear collective flow data in relativistic heavy-ion 
reactions [38,39,24,41] for densities up to five times saturation density. Also, 
the strength of the momentum-dependent isovector potential at po estimated 
from {Uneutron — U proton) /'^^ agrccs Very well with the Lane potential extracted 
from nucleon-nucleus scatterings and (p,n) charge exchange reactions with 
beam energies up to about 100 MeV [36,37,42,43]. 

We focus on the properties of spherically-symmetric, non-rotating, non-magne- 
tized neutron stars at zero temperature by solving the Tolman-Oppenheimer- 
Volkov equation. For the equation of state below about 0.07 fm~^, we use the 
results from Refs. [44,45]. Also, we assume that the neutron star consists of 
npept, matter, but does not contain any exotic components, such as hyperons, 
quarks, or Bose condensates. 

In Fig. 1 we display some of the basic properties of the EOS and the cor- 
responding neutron stars. Shown in the lower panel is the symmetry energy 
for X = 0, — 1 and —2, respectively. With x = the symmetry energy agrees 
very well with the prediction from Akmal, et. al. (APR) [46] up to about 5po. 
Around po, the x = EOS can be well approximated by Ef^{p) ^ 32(p/ po)°'^- 
With X = —1, the E^~^^{p) ^ 32(p/po)^'^ is closer to predictions of typical 
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Fig. 1. Neutron star mass as a function of the central density, and the proton fraction 
for beta-equilibrated matter and symmetry energy as a function of density for the 
EOS with a; = 0, —1, and —2. The dotted lines give the corresponding results for 
the APR EOS. 



relativistic mean field models [17]. 



The middle panel shows the proton fraction, function of density, and 

the top panel gives the mass of a neutron star as a function of the central 
density. For Xp below 0.14 [47], the direct URCA process does not proceed 
because energy and momentum conservation cannot be fulfilled. The proton 
fraction is sensitive to the slope of the symmetry energy [9,10,11,12]. For the 
X = —1 and X = —2 EOSs, the condition for direct URCA is fulfilled for nearly 
all neutron stars above 1 Mq. For the a; = EOS, the minimum density for 
direct URCA is indicated by the vertical dotted line, and the corresponding 
minimum neutron star mass is indicated by the horizontal dotted line. For 
the X = EOS, neutron stars with masses above 1.39 MqwIU have a central 
density above the threshold for the direct URCA process. 

Isospin diffusion in heavy-ion reactions is the re-distribution process of isospin 
asymmetries carried originally by the colliding partners. The degree and rate of 
this process depends on the relative pressures of neutrons and protons, namely 
the slope of the i?sym(p)- It is harder for neutrons and protons to mix up with 
a stiffer Esy^ip), leading to a smaller /slower isospin diffusion. Moreover, the 
distribution of the isospin asymmetry versus density during heavy-ion reac- 
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tions is completely determined by the Esym{p)- As an illustration, shown in 
the insert of Fig. 2 is a snapshot at 20 fm/c of the correlation between the lo- 
cal isospin asymmetry and density in the ^^^Sn+^^^Sn reaction using the two 
extreme density dependences of the £'sym(p). In this work, our calculations 
of nuclear reactions are performed using the latest isospin and momentum- 
dependent transport model using in-medium nucleon-nucleon cross sections 
consistent with the corresponding single particle potential [30] . With the very 
stiff symmetry energy of a; = —2, a very neutron-rich dilute cloud surrounds 
a more symmetric denser region up to 1.6po- With the very soft symmetry 
energy of x = 1 (it first rises then starts decreasing with the increasing p 
above about 1.3po [35], mimicking one of the results in Ref. [48]), however, the 
isospin asymmetries at both very low and very high densities are higher than 
the average asymmetry of the reaction system. The observed inverse relation- 
ship between the 6{p) and Esymip) is consistent with the well-known isospin 
fractionation phenomenon first predicted based on the thermodynamics of 
asymmetric matter [49,50]. Similar to neutron skins in heavy nuclei, neutron- 
rich clouds are dynamically generated in heavy-ion reactions via the isospin 
diffusion. This indicates that the same underlying physics is at work [29]. 

In the following, we examine the strength of the isospin diffusion and the 
thickness of neutron skin in ^osp]-, ^ function of the slope parameter L = 
3po{dEsym/dp)po. The degree of isospin diffusion in the reaction of A+B is 
experimentally measured by using [51] 



20f +B - 0f+^ - Of +^ 
0,^+^ - Of 



where 0/ is any isospin-sensitive observable. The frequently used ones in- 
clude the neutron/proton ratio of pre-equilibrium nucleons, ratios of light 
mirror nuclei and the isospin asymmetry of projectile-like fragments. They 
all give essentially the same result [27,52]. By construction, the value of Ri 
is 1 (—1) for the symmetric A + A [B + B) reaction. If a complete isospin 
equilibrium is reached in the asymmetric reaction A + B a.^ a. result of isospin 
diffusion the value of Ri is about zero. The Ri also has the advantage of re- 
ducing significantly its sensitivity to the symmetric part of the EOS. Shown 
in Fig. 2 are the strength of isospin diffusion 1 — Ri calculated using the trans- 
port model[30] and the size of neutron skin dRnp in ^osp]^ calculated using 
the Skyrme Hartree-Fock with interaction parameters adjusted such that the 
same EOS is obtained [29]. The strength of isospin diffusion 1 — i?j is seen 
to decrease, while dRnp increases with the increasing L as one expects. The 
NSCL/MSU data 1 - Ri = 0.525 ± 0.05 implies that the L{X) parameter is 
constrained between 62.1 MeV (x=0) and 107.4 MeV (x=-l). This is consis- 
tent with the meaurement dRnp = 0.2 ± 0.04 fm [53] and also with several 
recent calculations [17,32,33,34]. However, presently available measurements 
of the neutron skin thickness using hadronic probes have large systematic 
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Fig. 2. The strength of isospin diffusion in the ^^^Sn+^^^Sn reaction and the size of 
neutron skin in ^ospb as a function of the slope of the symmetry energy, respectively. 
The insert is the correlation of the isospin asymmetry and density at the instant of 
20 fm/c in the reaction considered. 

uncertainties associated with the strong interaction. 

The corresponding mass vs. radius curves for these EOSs, as well as for APR 
(using the AV18+6v+\JIX* interaction) are given in Fig. 3. In addition the 
constraints of causality, the mass-radius relation from estimates of the crustal 
fraction of the moment of inertia (A/// = 0.014) in the Vela pulsar [54], 
and the mass-radius relation from the redshift measurement from Ref. [56] 
are given. Any equation of state should be to the right of the causality line 
and the AI/I line and should cross the z = 0.35 line. The horizontal bar 
indicates the inferred limits on the radius and the radiation radius (the value 
of the radius which is observed by an observer at infinity) defined as Roo = 
R/^l - 2GM/Rc^ for a 1.4 M© neutron star. 

Since all three calculations with x = 0, — 1 and x = —2 have the same com- 
pressibility {Kq = 211 MeV) but rather different radii, it is clear that the 
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Fig. 3. The mass-radius curves for x = 0, —1, and —2 and the APR EOS. The hmit 
from causality, the Vela pulsar, and the redshift of EXO0748 are all indicated. The 
inferred radius of a 1.4 solar mass neutron star and the inferred value of Roo are 
given. 

radius is indeed rather sensitive to the symmetry energy while the maximum 
mass is only slightly modified [9,10,11,12,55]. The APR EOS has a compress- 
ibility of Kq = 269 MeV but almost the same symmetry energy as with x = 0. 
We note that the APR EOS leads to a 16% higher maximum mass (1.9Mq to 
2.2M0) but only a 5% decrease in radius (12.0 km to 11.5 km) as compared 
to the results with x = 0. 

Since only EOSs with symmetry energies between a; = and x = —1 are con- 
sistent with the isospin diffusion data and measurements of the skin thickness 
of lead, we take them as representative of the possible variation in neutron 
star structure that is consistent with terrestrial data. The APR and the x = 
EOS have nearly identical symmetry energies and slightly different radii. Neu- 
tron star radii are strong functions of the symmetry energy but also contain 
contributions from the isospin-symmetric part of the EOS, especially at higher 
densities. Even though the compressibility of the APR EOS is larger than that 
of the X = EOS, the pressure is typically lower in the APR EOS at densi- 
ties just above saturation, giving the APR EOS a smaller radius by about 5%. 
Thus we take this 5% difference as representative of the remaining uncertainty 
in the symmetric part of the EOS and extend the minimum radius to 11.5 km. 
Neutron stars with radii larger than 13.6 km are difficult to make without a 
larger symmetry energy or compressibility [17]. We conclude that only radii 
between 11.5 and 13.6 km (or radiation radii between 14.4 and 16.3 km) are 
consistent with the x = and x = —1 EOSs, and thus consistent with the 
laboratory data. It is interesting to note that a radius of R=12.66 km was 
recently predicted for canonical neutron stars using a new effective interac- 
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tion calibrated by reproducing several collective modes of ^°Zr and ^'^^Pb [34]. 
This radius falls right in the range of our constraints. Our constraints on the 
radius are also consistent with the range of radii from the extensive analy- 
sis in Ref. [17] with only a few exceptions. The field-theoretical models from 
this reference which are outside our suggested range either have a relatively 
large symmetry energy at saturation density (> 36 MeV), or have very soft 
symmetry energies created by extremely strong non-linear couplings which are 
atypical of most relativistic mean field models. 

Our results suggest that the direct URCA processes is likely for stars with 
masses larger than 1.39 Mq, which is the limit obtained from the x = EOS 
in Fig. 1. This constraint nearly matches the constraint for the direct URCA 
process of 1.30 Mq obtained in Ref. [34]. This is markedly different, however, 
from the result from APR, which gives a large threshold for the direct URCA 
process (even though the symmetry energy is very similar to our a; = EOS). 

Does this constraint agree with present neutron star radius observations? The 
answer to this question is "yes". Assuming a mass of 1.4 Mq, the inferred 
radiation radius, -Roo, (in km) is 13.5 ±2.1 [3,4] or 13.6 ±0.3 [5] for the neutron 
star in u Cen, 12.8 ±0.4 in M13 [6], U.btH for XT in 47 Tuc [7] and 14.5t|^ 
in M28 [8], respectively. Except the neutron star in M13 that has a slightly 
smaller radius, all others fall into our constraints of 14.4 km < Roo < 16.3 km 
within the observational error bars that are often larger than the range we 
gave. 

While the Vela AI/I upper limit does not provide any new information, the 
fact that the z = 0.35 line does not cross our range of radii implies a mass 
larger than 1.4 Mq for EXO-0748 (the minimum mass would be about 1.7 
Mq corresponding to the dot in Fig. 3). This is larger than the canonical 1.4 
Mq neutron star mass, but is not unreasonable since this object is accret- 
ing [56] ^ . 

While estimates of radii based on astrophysical observations are still very chal- 
lenging, it is useful to compare our results with recent Chandra/XMM-Newton 
observations. Together with more refined observations, future heavy-ion ex- 
periments (some recent progress in Ref. [58]) with advanced radioactive beam 
facilities [59] and measurements of parity violating electron-nucleus scatter- 
ing [60] will allow us to pin down more precisely the EOS of neutron rich 
matter. This would allow tighter contraints on neutron star radii. On the 
other hand, a neutron star radius measurement outside of our prediction may 
indicate non-standard physics. 



In fact, while revising this work, a measurement of the mass 2.10 ± 0.28 Mq and 
radius 13.8 ± 1.8 km for this object was reported in Ref. [57], confirming that this 
object is likely more massive than 1.4 Mq. 
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